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We show that the interaction between vortices and sound waves in atomic Bose-Einstein con- 
densates can be elucidated in a double-well trap: with one vortex in each well, the sound emitted 
by each precessing vortex can be driven into the opposing vortex (if of the same polarity). This 
cross-talk leads to a periodic exchange of energy between the vortices which is long-range and highly 
efficient. The increase in vortex energy (obtained by numerical simulations of the Gross-Pitaevskii 
equation) is significant and experimentally observable as a migration of the vortex to higher density 
over just a few precession periods. Similar effects can be controllably engineered by introducing a 
precessing localised obstacle into one well as an artificial generator of sound, thereby demonstrating 
the parametric driving of energy into a vortex. 

PACS numbers: 03.75.Kk, 03.75.Lm 



I. INTRODUCTION 

In a quantum fluid, such as an atomic Bose-Einstein 
condensate (BEC) or superfluid Helium, vortices possess 
quantized circulation, synonymous with them being a 
topological defect in the macroscopic condensate phase. 
Quantized vortex lines, rings, lattices and tangles have 
been the subject of experimental study in superfluid He- 
lium for over 50 years [1 , and in which recent empha- 
sis has been on their role in quantum turbulence [2j [3] . 
Meanwhile, since the late 1990's, there has been fast- 
growing interest in vortices in Bose-Einstein condensates 
[U [5] , where the controllability of these gases has led to 
the generation of single vortices [6j [7] , giant vortices [8] , 
vortex dipoles [9] , soliton- vortex hybrids [10] and turbu- 
lent vortex tangles [11]. It is worthy of note that recent 
breakthroughs in imaging of both Helium [12 and BEC 
[T3] systems now enable the dynamics of quantized vor- 
tices to be monitored in real-time. 

The full nature of the vortex-sound interaction in quan- 
tum fluids is unclear [14 . The superfluid topology con- 
strains quantized vortices to disappear by annihilating 
with an oppositely charged vortex or vanishing at the 
edge of the system (where they annihilate with their im- 
age). In the limit of zero temperature and for a uniform 
condensate, sound waves are the low- lying excitations of 
the system and provide the only energy sink for vortex de- 
cay. For example, at zero temperature the reconnection 
of vortex lines [15] and the acceleration of a vortex line 
segment both generate sound waves [T6J [17] , dissipating 
the vortical energy. In the latter case, the acceleration 
that drives sound emission may arise from the influence 
of other vortices [TSJ [19] , Kelvin waves excitations of vor- 
tex lines [20] or the Magnus force in an inhomogeneous 
ambient density, e.g. in a trapped condensate [21 . The 
experimentally observed decay of vorticity at very low 
temperature in superfluid He [22j [23] is thought to be 
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primarily due to the Kelvin waves dissipation route, with 
reconnections playing a secondary role [3]. 

Less well understood is the inverse process, i.e. the 
absorption of sound by a vortex. Insight may be gleaned 
from vortex-sound interactions in fluid dynamics [14]. 
For example, an acoustic ray model has predicted that 
certain trajectories of sound wave can spiral into the vor- 
tex core, transferring energy to the vortical flow [24| [25]. 
While sound waves can induce the nucleation of vortices 
through the collapse of cavitating bubbles [26, 27], sound 
absorption by pre-existing vortices is not thought to play 
a significant role in homogeneous superfluid Helium sys- 
tems. However, sound absorption may become consider- 
able in atomic BECs due to their confined geometry. In- 
deed, the lack of sound-induced decay of a vortex precess- 
ing in a harmonically-trapped BEC has been attributed 
to reabsorption of the emitted sound [21 (although re- 
lated works [4j [16] predict that the sound emission may 
be prohibited due to the sound wavelength exceeding the 
system size). The harmonic nature of the trap appears 
key to supporting a sound-vortex equilibrium, with trap 
anharmonicities apparently leading to net vortex decay 
[2T| [28] . in analogy to dark solitons [29j [30]. However, 
it is difficult to resolve the interaction of sound and vor- 
tices in single trapped condensates [19] [21] due to their 
co-habitation in the trap. 

Analogous questions exist over the interaction of 
matter- wave dark solitons with sound. Like vortices, they 
radiate sound waves under acceleration J29[ [31] , become 
stabilised in harmonic traps [21] and can be parametri- 
cally driven by artificially generated sound waves [T8] . 
Recently, a double-well trap was shown to offer benefi- 
cial insight into the soliton-sound interaction [32] . With 
one soliton in each well, large-scale exchanges in energy 
between the solitons clearly demonstrated a long-range 
sound-mediated interaction. 

Here we will exploit a double well trap to study the 
emission and absorption of sound by vortices. We will 
begin by considering an idealized double trap geometry. 
The trap system and theoretical model are outlined in 
Section II. In Section III we consider how a single vortex 
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behaves in this system (while not the main results of our 
work, this an essential prerequisite to understanding the 
dynamics in later sections). In Section IV we progress to 
consider how two vortices, one in each well, interact via 
the exchange of sound waves or "cross-talk" . In Section V 
we replace one of the vortices with a moving obstacle, and 
explore how the sound generated interacts with the re- 
maining vortex. In Section VI we discuss our theoretical 
findings. In Section VII we demonstrate the same quali- 
tative behaviour in an experimentally-achievable double 
trap geometry. Finally, in Section VIII we draw conclu- 
sions of our work. 



II. THEORETICAL FRAMEWORK AND TRAP 
SET-UP 

We consider a BEC at ultracold temperatures such 
that thermal and quantum fluctuations can be neglected 
and that the system is well parameterised by a mean- 
field order parameter \I/(r, t) which satisfies the three- 
dimensional Gross-Pitaevskii equation (GPE) [55] . We 
assume a quasi-2D geometry, whereby harmonic trapping 
is sufficiently tight in one dimension, taken here to be the 
^-direction, to freeze out the corresponding dynamics. 
Then, using the decomposition \I/(r,£) = i/j(x,y,t)i/j z (z), 
one can integrate out the time-independent axial com- 
ponent i\) z {z) from the 3D GPE. The transverse order 
parameter ip(x,y,t) then satisfies the 2D GPE, 

ihd t ip = (-^v 2 + V(x,y) + <#| 2 - /i) ^ (1) 

where V(x, y) is the axial trapping potential and m is the 
atomic mass. The 2D chemical potential fi is related to 
the 3D chemical potential // via \i = [i f — huj z /2, where uo z 
is the trap frequency in the axial direction, s- wave atomic 
interactions, of length a, give rise to the nonline ar term 
with coefficient g = 2\/2iiTi 2 a/ ml z , where l z = ^J?i/muj z 
is the harmonic oscillator length of the frozen dimension. 

The complex order parameter ip(x, y, t) can be written 
as i/j(x,y,t) = n(x, y, t) exp[i0(x, y, t)] , where n(x,y,t) 
and (j)(x,y,t) are the distributions of atomic density and 
phase, respectively. Furthermore, the phase defines the 
fluid velocity v = (Ti/m)V(j). In 2D vortices are singular 
points about which the phase wraps around by an integer 
multiple of 2tt and the condensate flows azimuthally. The 
density is pinned to zero at the central point creating a 
well-defined vortex core which relaxes to its unperturbed 
value at a distance of the order of the healing length 
f Tt/y/mng. 

We will initially consider an idealized double well sys- 
tem consisting of two connected harmonic traps, 
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FIG. 1: Vortex-free density profile at y = along the x-axis 
(black solid line) in the double harmonic potential (grey solid 
line) with a barrier height at x = of Vb = 0.6/x . The 
corresponding Thomas- Fermi profile utf = {l^ — V)/g is also 
shown (black dashed line). 



[T] The barrier separating the trap, which has a mini- 
mum height of Vb = \muJ 2 x 2 c , determines the connectiv- 
ity between the wells and thus the degree to which sound 
and vortices can propagate between wells. The transfer 
of sound waves, which have energy of order /i, between 
the wells will be possible for Vq < /J> and prohibited for 
Vb ^> /i. The capacity for vortices to propagate between 
wells will depend additionally on the energy of the vortex. 

While this trap is not directly achievable in experi- 
ments due to its sharp feature at x = 0, it is convenient 
to consider theoretically: it identifies the main physical 
effect in a "clean" manner and allows us to draw on the 
established knowledge for vortices in harmonic traps. We 
will later demonstrate that the same dynamics persist in 
experimentally-achievable set-ups. 

The 2D GPE is solved numerically using the Crank- 
Nicholson method [34]. Over a typical simulation, the 
relative change in norm and energy AN/N and AE/E 
are of order 10 -6 , i.e. the solution is numerically well 
converged. The vortex-free ground state, with density 
profile nyp(x^y)^ is found by propagating the 2D GPE 
in imaginary time. Vortical states are imposed by forcing 
the phase distribution (during imaginary time propaga- 
tion) to, 



,y) = JJ^arctan 



y-yj 



(3) 



Each trap is circularly symmetric with frequency uo and 
displaced from the origin by d=x c , as illustrated in Fig. 



where i is the index of a vortex with charge qi located 
at (xi,yi). Since multiply-charged vortices are energet- 
ically unfavourable compared to multiple singly-charged 
vortices [TJ [55] , we shall only consider singly-charged vor- 
tices = 1 (relative polarity may change). We will infer 
the change in energy of vortices in our system through 
changes in their position, e.g. a drift to lower density 
signifies a decrease in the vortex energy. This is exactly 
what would be done in reality, since the vortex energy 
cannot be directly measured experimentally. 

We assume units in which length, speed and energy 
are expressed in terms of the 2D healing length £ = 
Tij yjmn^g, speed of sound c = y^u^g/m and chemical 
potential /i = nog, where no is the peak density. 
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The ratio ji/ftuo specifies the nature of the conden- 
sate. For ji/fiuj <C 1 the trap dominates and the ground 
state will approximate the gaussian harmonic oscillator. 
For h/Tiuj ^> 1 the interactions dominate and lead to a 
broad condensate profile. Then the kinetic energy of the 
ground state, which depends on the gradient of the den- 
sity, becomes negligibly small. Under this Thomas-Fermi 
(TF) approximation, the density takes the analytic form 
n TF = ( f i-V)/g [33]. 

We focus on a system with ji/ftuj = 10. An example 
density profile is shown in Fig. [I] It is closely matched 
by the TF prediction, with the only significant devia- 
tion arising at the condensate perimeter and the bar- 
rier. There the density gradient is not negligible, result- 
ing in a smoothening of the density over a lengthscale 
~ The TF approximation predicts a condensate ra- 
dius Rtf = — 14.14£. While we hereafter 
express length in terms of healing length, this can be 
trivially related to the trap harmonic oscillator length 
/ho = y^h/muj via Zho = >/l0£- 



III. SINGLE VORTEX 

We first explore the dynamics of a single vortex within 
the double trap system (J2|. A vortex tends to follow 
a path of equipotential through a trapped condensate, 
e.g., precessing around the trap centre in a single har- 
monic trap [55] . In the double trap we can addition- 
ally anticipate a regime in which the vortex can follow 
a dumbbell-shaped path around both wells. As is well 
known for single harmonic traps, the energy and angular 
momentum associated with the vortex increases as the 
vortex is moved to higher density, i.e. towards the centre 
of the well [1 [55]. 

We place a vortex in the right-hand well at a position 
(x c ,y c + r v ), where (x c ,y c ) = ([2V /muj 2 } 1/2 , 0) is the 
origin of the right-hand trap and r v is the initial offset 
of the vortex from the well centre. NB the dynamics are 
insensitive to the direction of the off-set. The ensuing 
dynamics are sensitive to the size of the vortex displace- 
ment r v and the inter-trap barrier Vq. We numerically 
evolve the vortex dynamics over a long simulation time 
(5000 (£/c)) within this parameter space. Note when the 
vortex is initially positioned very close to BEC edge (typ- 
ically within one healing length of .Rtf) its evolution is 
indistinguishable from surface excitations that are gen- 
erated. This is a general consequence of being initiated 
in the low density periphery and occurs even in isolated 
harmonic traps. As such we do not present the vortex 
dynamics at such extreme positions. 

The phase diagram for the vortex dynamics is plotted 
in Fig. [2] separating regions of 'stable' vortex dynamics 
(dots) from those where there is no vortex within the sys- 
tem in the 'final' simulated state (crosses), as discussed in 
detail below. Note that we observe a qualitatively similar 
phase diagram for a larger (more TF-like) condensate. 



A. Stable dynamics 

If dissipationless, we would expect the vortex to always 
remain in the system. In Fig. [2] we see regimes where this 
is true (dots) but also where the vortex is unstable and 
ultimately leaves the system (crosses). Stable vortex mo- 
tion is promoted for large Vo, since the well then behaves 
likes an isolated harmonic trap, and for low vortex radii, 
since the vortex does not feel a strong effect from the far 
trap. Under this stable motion the vortex dynamics is 
akin to that in a single harmonic trap [36 : it precesses 
around the well centre with approximately constant ra- 
dius and generates a collective motion of the background 
condensate of low amplitude (~ 5%no). In Fig. [3] we 
show a typical snapshot of the condensate. The den- 
sity distribution consists of two weakly-connected circu- 
lar condensates and the vortex appears as a hole (white 
spot) in the right-hand well. By subtracting the time- 
independent vortex-free density nyF from this density 
profile, the collective excitation become clearly visible 
(Fig. [3]). The precession frequency of the vortex is ~ 0.2a; 
for small displacements, and increases with r v , in good 
agreement with analytic predictions for vortex precession 
in a single harmonic trap in the TF regime [57] . 
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FIG. 2: Phase diagram in r v — Vb space for the single vortex 
case showing whether the final state of the system (after a 
long simulation time of 5 x 10 3 (£/c)) is a single vortex (dots) 
or the vortex- free state (crosses). The light-shaded region is 
the cross-over regime (case I), in which the vortex passes into 
the far well. The dark-shaded region is the inductive regime 
(case II), in which the initial vortex induces other vortices. 
The lines are the TF predictions for the onset of cross-over 
dynamics t v /Rtf — — n m in/no) using the TF predic- 

tion (solid line) and the numerical values (dashed line) for 
rimin- The observed deviation between these lines and the 
light/dark grey boundary acts as an indication of the impor- 
tance of sound emission in the actual vortex dynamics. The 
TF radius of each well is Rtf = 14.14£. 
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FIG. 3: (Left) Density n(x,y) for the single vortex scenario 
with y v = 1£ and Vb = 0.9/i at a time of 1000 £/c. Black 
(white) corresponds to peak (zero) density. (Right) Renor- 
malised density n(x,y) — 71vf(^, y) for the same data as above. 
The colorscale is ±10%no- Each box is of size 64£ x 32£. 



B. Unstable dynamics 



In cases where the vortex eventually decays from the 
system (crosses in Fig. |2| its initial dynamics falls into 
one of two cases. In case I, the dynamics is charac- 
terised by the vortex crossing over into the adjacent well, 
whereas in case II it is characterised by inducing a mirror 
vortex in the adjacent well. These effects most commonly 
become manifested in the first precession of the vortex in 
the trap, but in a minority of cases they may arise after 
several precessions. 



1. Case I: vortex crossover 



Case I (vortex crossover) arises in the lightly shaded 
region in Fig. [2] The vortex can be expected to travel 
between the wells when the local potential of the vortex 
(the value of the potential at the vortex core) exceeds 
the inter-well barrier. This will tend to occur for large 
vortex off-sets and a weak inter-trap barrier. Put quan- 
titatively, one would expect cross-over to occur when the 
vortex density depth n v is less than or equal to the min- 
imum density at the barrier n m i n (one can picture this 
as when the vortex can just squeeze through the bar- 
rier). If dissipationless, the density depth of the vortex 
will retain its initial value, which we can approximate via 
the TF prediction n v = no(l — rl/R 2 ). This is a robust 
approximation provided that the vortex position is away 
from the edge and the barrier, for which the TF density 
agrees with the actual density profile to within O.Olno- 
This gives the criteria r v /R > ^/l — n m i n /no for cross- 
over to be possible. We can first approximate n m i n via 
its TF prediction nj^f n = no(l — Vb//i), giving the solid 
line. However, the onset of cross-over dynamics occurs 
at considerably lower r v . We can expect some deviation 
to arise from the inaccuracy of using the TF approxima- 
tion for n m i n : the TF approximation underestimates the 
density at the point of the barrier, as evident in Fig. [I] 




FIG. 4: Crossover regime dynamics: Snapshots of (left) 
density n(x,y) and (right) renormalized density n(x,y) — 
nvF(x,y) within the crossover regime, for r v = 8£ and 
Vb = 0.6/x, at various times. (Spatial and color scales are 
the same as in Fig.pl. 



If we instead use the actual value of the ground-state 
density at the barrier, we obtain the threshold shown 
by the dashed line in Fig. [2] This lowers the prediction 
for r v , but only slightly and the prediction still remains 
considerably greater than the observed threshold. This 
anomaly is likely to arise from the radiation sound from 
the vortex. This will have two implications for promot- 
ing cross-over dynamics. Firstly the dissipating vortex 
will move to lower densities/greater radial position. Sec- 
ondly, the emitted sound generates collective modes of 
the BEC which will cause the density at the barrier to 
become time- dependent and may promote vortex cross- 
over during "high tide" . 

In this crossover regime, the ultimate fate of the vor- 
tex is to decay. In Fig. [4] we present an example. The 
precessing vortex approaches the barrier (t = 30(£/c)) 
and upon traversing it (t = 60(£/c)), decays into a high- 
amplitude curved pulse of sound (t = 67.5 (£/c)). The 
sound pulse reflects off the far left-side of the trap and 
propagates back through the trap (t = 120 (£/c)). Fol- 
lowing many reflections and diffractions in the trap, the 
sound pulse becomes randomised, ultimately forming an 
isotropic sound field (t = 375 (£/c)). 

In other cases, the vortex undergoes a more gradual 
decay, passing many times between the wells. NB more 
generally, we can observe some cases where the vortex 
undergoes stable cross-over dynamics, as we will see later 
in Section VII. 
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FIG. 5: Inductive regime dynamics: (a) Density n(x,y) (left 
column) and renormalised density n(x,y) — 71vf(^,2/) (right 
column) at various times. Parameters: r v = 5£ and Vb = 0.9/i 
(box size and colorscale are the same as in Fig. [3}. (b) Evo- 
lution of the radial position of the vortex r v ; arrows indicate 
the points at which vortex induction takes place. 



vortex in the right-hand well. As a result there can arise 
periods of time where multiple vortices can appear in the 
system. For all cases we observe the growth of a tempes- 
tuous sound field during the vortex motion, set up by the 
sound emission from the accelerating vortices. All of the 
vortices ultimately dissipate and disappear into an ener- 
getic and isotropic sound field (Fig.[5jb) at t = 1540^/c). 

Up until the point when it disappears, the original vor- 
tex drifts outwards, as shown in Fig. j5^b) . Most of the 
noise in r v (t) is due to the buffetting effect that the sound 
field has on the vortex. However, the sizeable jumps in 
r v at t ~ 300 and 700 (£/c) (highlighted by arrows) oc- 
cur as a vortex is induced in the far well, indicating the 
transfer of energy from the original vortex to create the 
new one. 

This region of the parameter space occurs for a range 
of barrier heights in the vicinity of Vb = /i, for which 
the condensate channel is of low density. The energy to 
create a vortex depends on the local density. As the bar- 
rier height is reduced below /i, the density in the barrier 
region increases and we have confirmed with numerical 
simulations that the energy cost to create a vortex here 
increases sharply. 



IV. CROSS-TALK OF TWO VORTICES 

We now extend to the case where there is initially a 
vortex in each well to examine the possibility of sound- 
induced interaction between vortices. We continue to 
employ the idealized double harmonic trap geometry; 
we will demonstrate the same phenomena in a realizable 



double trap geometry in Section |VII| The additional vor- 
tex (denoted "vortex 2") is created at the centre of the 
left-hand trap. Since it feels the velocity field of vortex 
1 it will is not perfectly stationary. The displacement 



2. Case II: vortex induction 

The second case of unstable dynamics (case II: vor- 
tex induction), characterised by the initial vortex induc- 
ing a second vortex in the far-well, arises in the dark 
shaded region in Fig. [2] Snapshots of a typical evolu- 
tion are shown in Fig.j5] As the vortex precesses close 
to the adjoining well (t = 300£/c), a mirror vortex be- 
comes excited on the opposite side of the inter-well bar- 
rier. This vortex has the opposite charge to the origi- 
nal vortex and precesses in the opposite direction around 
its trap (t = 415£/c). This early induction is the char- 
acteristic of this regime of dynamics. The subsequent 
dynamics can vary widely. In the example, the initially- 
induced vortex disappears (t = 500^/c), a second vortex 
is induced (t = 720^/c) and then the original right-hand 
vortex disappears (t = 920^/c). In other cases the origi- 
nal vortex may disappear as it creates its mirror vortex, 
while sometimes the induced vortex may itself induce a 
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FIG. 6: Phase diae; ram m Tv — Vb of the final state (after 
5000 (£/c)) of the two vortex system for (a) same polarity and 
(b) opposite polarity of the vortices. The final state is either 
two stable vortices (dots), one stable vortex (pluses) or no 
vortices (crosses). The second vortex is initially positioned at 
the centre of the left-hand well. Light and dark shading indi- 
cates case I (crossover) and case II (induction) of the unstable 
dynamics. 
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FIG. 7: Vortex dynamics in the two- vortex system with (a) the same polarity and (b) different polarity. Parameters: n(t = 
0) = 4£, Vb = 0.6//. Black (green/dark grey) lines correspond to vortex 1 (2). (c) and (d) plot the corresponding vortex radii, 
and additionally include the evolution of a single vortex in an isolated harmonic trap (light blue/grey line) with initial radius 
ry (t = o) = 4£. 



of vortex 1, denoted r v , and the height of the inter- well 
barrier, Vb, are again key to the dynamics. The latter pa- 
rameter now additionally controls the transfer of sound 
between the wells, since sound waves have an energy of 
around fi. We can also expect sensitivity to the relative 
polarity of the vortices. Figure [6] shows the phase dia- 
gram of the system for (a) vortices of the same polarity 
and (b) vortices of opposite polarity. 

The phase diagrams are similar to the single vortex 
case (Fig. |5| and the vortex polarity only has a small 
effect on the final state of the system. There are two sta- 
ble regions, one for weak barriers and low displacements, 
and the other for large barriers. We also see a cross-over 
regime for weak barriers and large vortex displacements 
in which vortex 1 tends to cross into the other well, and 
an induction regime for barrier heights centered around 
\i in which vortex 1 induces another vortex in the oppo- 
site well. Where vortex instability does occur, the end 
state is either the persistence of a single vortex or no 
vortices at all. We cannot make any general comments 
about what favours these two end states: the dynamics 
that can develop are sufficiently complex that the end 
state is not readily deterministic. 

Note that in Fig. [6^b) there exists stable dynamics in 
the cross-over regime. In this case vortex 1 perpetually 
traverses both wells in a stable manner while vortex 2 
remains localised in the centre of its well. 

In order to investigate the cross-talk between vortices, 
that is their sound- mediated interactions, we focus on 
the stable regimes and in particular the stable area oc- 
curring for low barrier heights, for which one can expect 
unimpeded motion of sound between the wells. Figure [7] 



present the vortex dynamics for r v = 4£ and Vo = 0.6/i. 
When the vortices have the same charge (Fig. |7ja)) the 
vortices periodically spiral inwards and outwards, and 
remain out of phase with each other. The change in 
vortex radial position (Fig. [7^c)) is large (the vortices 
oscillate between the trap centre and a radius of around 
6£ = 0.43i?TF), demonstrating a significant transfer of 
energy between the vortices. The energy exchange oc- 
curs over a timescale of ~ 2000 (£/c), which is around 
8 precessions of the vortex in the trap (the precession 
period is ~ 270 (£/c)). In the analogous situation for 
the dark soliton, the energy transfer is much slower, oc- 
curring over many tens of soliton oscillations [32 . The 
periodic change in vortex radius (Fig. (7^c)) is composed 
of sharp steps, suggesting that the energy emission and 
absorption of the vortices occurs in packets rather than 
continuous. 

In Fig.[7jc) we also present, for comparison, the evolu- 
tion of a single vortex in an isolated harmonic trap (light 
blue/grey line). The vortex maintains an approximately 
constant radius, with only small-scale modulations due 
to its re-interaction with its emitted sound and collective 
modes [36] . 

When the vortices have opposing polarity (Fig. (7^b) 
and (d)), no significant transfer of energy is observed 
and the vortices precess with approximately constant ra- 
dius. The vortex motion does undergo modulations but 
these are not out-of-phase and so do not constitute a di- 
rect exchange of energy. Rather the modulations arise 
from the back-action of the randomised sound field on 
the vortices. These modulations are larger than the cor- 
responding ones experienced by a single vortex in a har- 
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monic trap, indicative of the greater sound density in the 
double well system. 

The significant transfer of energy between like-charged 
vortices and insignificant transfer between unlike-charged 
vortices is consistent across the stable region of the pa- 
rameter space at low barrier heights. In the region of 
stable dynamics at high Vq we do not observe a signif- 
icant transfer in energy between the vortices since the 
transfer of sound across the barrier becomes prohibited. 

It would appear that the vortices periodically drive en- 
ergy into each other via their emitted sound. However, 
we must rule out that this effect arises from the long- 
range interaction between vortices due to their superim- 
posed velocity fields. To this aim, we will next attempt 
to drive a vortex via sound waves generated by artificial 
means - a moving localised barrier. 



V. PRECESSING OBSTACLE 

We now replace the left-hand vortex with a precessing 
obstacle. Note that we continue to employ the idealized 
double harmonic trap; the same effects will be demon- 
strated in an experimentally realizable trap in Section 
The obstacle corresponds to a time-dependent po- 
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(4) 



{x + x c - x oh (t)} 2 + {y- y b(t)} 2 



tential, 

V Q h(x, y,t) = A oh 
x exp 

where, 



^ob(^) = ±^ob sin(cj b£ + </>ob) 
Vob(t) = r oh cos(u oh t + 0ob). 

This represents a Gaussian-shaped barrier, of amplitude 
A Q b and width a, moving in a circular path around the 
centre of the left-hand well with radius r Q b = + y^ h . 
Such an obstacle can be induced experimentally via a 
blue-detuned laser beam [9, 38 j. For the "-" ("+") case 
above, the barrier moves in the same (opposite) direction 
as the vortex. <p \^ represents the relative phase between 
the initial position of the obstacle and the vortex (e.g. for 
0ob = the obstacle and vortex start at mirror-opposite 
positions in their wells). 

We attempt to generate sound that mimics the sound 
generated by a vortex and will thus employ precession 
frequencies similar to those of a vortex (~ 0.23co> for our 
system). This is sufficiently low to prevent the obstacle 
exceeding the superfluid critical velocity and nucleating 
vortices [9, 38-40 . We have confirmed via numerical 
simulations of the precessing obstacle (with no vortices 
imposed) that its only effect is to generates spiral sound 
waves with typical amplitude ~ l%n$. 

We mimic the system parameters employed in Fig. [7] 
[r v = 4£ and Vo = 0.6/i] and additionally employ an 
obstacle potential defined by width <j = 2£, amplitude 




t/fi/c) 



x 10 



FIG. 8: Evolution of the vortex radius under the influence of 
a driving obstacle. Black line (B): the obstacle moves in the 
same direction as the vortex. Light grey line (A): the obstacle 
moves in the opposite direction to the vortex. Dark grey line 
(C): the barrier moves in the same direction as the vortex 
but its motion is terminated at t=8800(£/c). Parameters: 
n(t = 0) = 4f, a = 2f, r ob = 4f, ob = and cj ob = 0.234cj 
and Vo — 0.6/x. 



A b — phase o b = 0, path radius r Q b = 4£ and 
frequency cj b = 0.234cj (these parameters are chosen 
as they optimise the energy transfer, as shown below). 
First consider the case when the barrier moves in the 
same direction as the vortex (solid black line in Fig. [8|. 
We clearly see that the vortex radial position oscillates 
in time, i.e., there is periodic driving of energy into the 
vortex. This is analogous to the predicted parametric 
driving of a matter- wave dark soliton [18] . The energy 
transfer due to the obstacle is more gradual than that 
induced by another vortex (i.e., compare to Fig. [7]), with 
the exchange occurring over a time of ~ 12000 (£/c) or 
~ 40 precessions. Conversely, if the barrier moves in the 
opposite direction to the vortex (solid green/grey line) 
no driving is observed (although a small outward shift of 
the vortex suggests a small loss of energy in this case). 
This is qualitatively the same behaviour as we observed 
earlier when the obstacle in replaced by another vortex. 
However, the timescale of the energy transfer is consider- 
ably slower than with the obstacle. This shows that the 
vortex is a much more efficient source of sound to drive 
another vortex. 

Returning to the case where the obstacle moves in the 
same direction as the vortex, we find that if the obstacle 
motion is terminated (becoming a stationary obstacle) 
when the vortex is at its minimum radial position (dotted 
black line in Fig. [8|, it maintains this minimal radial 
position (subject to some oscillations about this radius). 
In this manner we can parametrically drive net energy 
into the vortex. Indeed, if we terminate the obstacle 
motion at any point the vortex approximately maintains 
that radial position. 

The minimum radius achieved r m i n provides a measure 
of the maximal amount of energy driven into the vortex. 
In Fig. [9] we show how this quantity varies as a function 
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FIG. 9: Minimum radial position achieved by the vortex un- 
der the influence of a precessing Gaussian barrier during a to- 
tal simulation time of 2 x 10 4 (£/c). The vortex is initially at 
r v = 4£ and the inter- well barrier height is Vb = 0.6/i. The ob- 
stacle parameters are based around the values cj G b = 0.234cj, 
A h = /i, cr = 2£, r Q b = 4£ and o b, and each is varied 
in turn from plot (a)-(e). In (a), we present results for co- 
rotating (solid line with dots) and anti-rotating (dashed line 
with crosses) system. In (b)-(e) we consider only the co- 
rotating system. In each plot the initial vortex radial position 
is shown (dashed horizontal line). In (a) we show the preces- 
sion frequency of the undriven vortex (grey vertical line). 

of the obstacle parameters lj b,V b, ^ob? cr and o b- For 
each case we observe a clear resonance that minimises 
r m i n and maximises the energy that can be driven into 
the vortex: 

• o; b [Fig. [9^ a)]: When the obstacle direction is the 
same as the vortex (black line/dots) we see a clear 
frequency resonance, with maximal energy being 
driven into the vortex for a precession frequency 
close to that of the undriven vortex (vertical dotted 
line). For the opposite case, the minimum radius 
achieved undergoes no significant resonance. For 
this reason, we subsequently only consider the co- 
rotating case. 

• Vob [Fig. [9^b)]: r m i n is minimised for an obstacle 
height of ~ fx. 

• a [Fig.^c)]: A barrier width a ~ 2£ best promotes 
energy driving. 

• r Q b [Fig.|9jd)]: The barrier position shows a double 
resonance at 4£ and 8£. 



• (/>ob [Fig. [9je)]: Maximal energy driving occurs 
when the vortex and obstacle begin in phase or with 
the obstacle slightly lagging the vortex. 

VI. DISCUSSION 

A. Cross-over and induction dynamics 

Our precursive study of a single vortex in a double 
trap revealed cross-over and inductive dynamics. The 
generation of a vortex state in a neighbouring well of 
a double- well trap has been studied in [41]. However, 
there the vortex-containing BEC partially tunnelled into 
the adjacent empty well, and so is a distinct creation 
mechanism to the one observed here. While these regimes 
are interesting in their own right, from the perspective 
of exploring cross-talk and parametric driving of vortex 
they pose limitations on the parameter space in which 
vortices are stable. 



B. Cross-talk between vortices 

Our main results demonstrate cross-talking between 
two vortices. When the vortices have the same polarity 
and different initial positions, they undergo periodic ex- 
changes of energy with each other, transferred via sound 
waves. We may view this as the driving of a vortex by 
sound from another vortex. The transfer of energy is 
rapid: a full energy cycle typically occurs within 10 vor- 
tex precessions. 

Sound is emitted from each vortex due to its accelera- 
tion in the trap. This emitted sound has a quadrupolar 
radiation pattern and forms an outwardly propagating 
spiral wave [21] that carries angular momentum of the 
same orientation as the originating vortex (to conserve 
angular momentum). Upon passing into the opposing 
well it imparts some of its energy and angular momen- 
tum to the opposing vortex. If the angular momentum 
carried by the sound waves is of the same orientation as 
the receiving vortex (which is to say that the vortices are 
of the same polarity), it serves to increase the angular 
momentum and energy of this vortex, causing it to move 
towards the high density trap centre. By the same simple 
argument, one would expect that when the vortices are 
of opposite polarity, the energy and angular momentum 
of the opposing vortex would be reduced. However, no 
exchange of energy is observed when the vortices have 
opposing polarity. We will return to this anomaly below. 
The cross-talk between vortices is analogous to that oc- 
curring for dark solitons in the ID analog of this system 
[32] . However, where the energy transfer for solitons oc- 
curs over many tens of oscillations, the energy transfer 
for vortices can occur in less than 10 oscillations. Thus 
the vortex system appears a beneficial platform to exper- 
imentally observe and explore these related sound phe- 
nomena. 
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C. Parametric driving of a vortex 

By moving a localised Gaussian obstacle through one 
well, we are able to generate angular- momentum carry- 
ing sound waves. When this angular momentum is of 
the same sign as the vortex, this acts to drive energy 
into the vortex, but when it of opposite sign it is essen- 
tially invisible to the vortex, as seen above for the vortex- 
vortex transfer. Under constant driving, the vortex en- 
ergy oscillates periodically and the oscillation amplitude 
can be resonantly tuned via the obstacle parameters. The 
driving is most effective when the time-dependent den- 
sity perturbation created by the obstacle closely matches 
that of the vortex, i.e. the same width, precession fre- 
quency and amplitude. However, the vortex-obstacle en- 
ergy transfer is considerably slower than vortex-vortex 
transfer, suggesting that the most "natural" sound to 
drive a vortex is that from another vortex. 

An analogous set-up was studied in [42], except the 
vortex and obstacle resided in the same trap. There the 
precessing obstacle periodically formed a vortex at the 
condensate edge and drove it into the centre. The dy- 
namics were interpreted as nonlinear Rabi cycling from 
the ground state to the first vortex state. Similarly, we 
may interpret our observations as oscillations between 
a high and low energy vortex state. Our results suggest 
that sound waves may have played a central role in trans- 
ferring energy from the obstacle to the vortex in [42] . 

Our findings mirror those relating to driving of dark 
matter- wave solitons in [18]. There, ID oscillating pad- 
dles were employed to impart linear momentum to the 
condensate sound field which in turn led to energy being 
driven into the dark soliton. Resonant driving occurred 
when the paddles were oscillated at close to the frequency 
of the undriven soliton motion. Furthermore, when the 
drive was switched off the soliton maintained that energy 
(in the absence of phenomenological dissipation), subject 
to oscillations from the background density excitations. 



D. Insight into sound absorption in trapped and 
homogeneous systems 

Using an acoustic ray model, Nazarenko et al. [24j [25] 
find that certain trajectories of sound waves incident on 
a vortex line can spiral into the vortex core, imparting 
their energy to the vortical flow. Despite this and our ob- 
servations herein, sound absorption is not cited to play a 
significant role in superfluid Helium systems, e.g. quan- 
tum turbulence [2j [3] . 

Consider an infinite homogeneous system containing a 
tangle of vortex lines. Within this, consider an element of 
vortex line emitting a pulse of sound due to some acceler- 
ation. From the point of emission, the sound spreads out 
radially. This rapidly dilutes its energy and momentum 
density, and thus its capacity to influence a line element 
in its path. Superimpose many such events from many 
randomly oriented vortex lines and the combined sound 



field will be tends towards being isotropic with no net an- 
gular momentum. As such it is not surprising that sound 
absorption is insignificant in these systems. 



Now consider our trapped BEC. Sound radiated out- 
wards by an accelerating vortex will eventually reflects off 
the trap wall and become partially focussed towards the 
trap centre. While the effectiveness of the focussing will 
vary with vortex position and trap shape, it will nonethe- 
less lead to a greater and more sustained sonic energy and 
momentum density than in a homogeneous system, which 
a suitably placed vortex may be able to gain from. The 
focussing effect will, however, be shortlived as the sonic 
angular momentum will become randomised after several 
reflections in the trap. We see this in Fig. [4] where a large 
sound pulse with net momentum is rapidly randomised 
into an isotropic sound field. This further explains the 
rapid equilibration of the vortex energy when the driv- 
ing is terminated. These simple arguments suggest that 
sound absorption may play a significant role in trapped 
condensates. 



Irrespective of the source of the sound, we observe 
transfer of energy between sound and vortex only when 
they have the same sign of angular momentum. We 
may interpret this in terms of the acoustic ray picture 
[24] [25] which predicts that trajectories of sound spiral 
into the vortex core and dump their energy and momen- 
tum. There the sound was modelled as plane waves while 
in the BEC the incident sound is a spiral wave carrying 
angular momentum. The crucial aspect is that, when 
the sound waves have the same angular momentum as 
the vortex, they naturally wrap around the vortex, pro- 
moting them to spiral into the core. Conversely, when 
the sound has opposite angular momentum, it will tend 
to be deflected and repelled by the vortex, unable to im- 
part its angular momentum. In essence, the vortex is 
invisible to sound waves of opposite angular momentum. 



In our quasi-2D system, the vortex line is rectilinear 
and can only raise its energy by moving to regions of 
higher density. In 3D system, the vortex can addition- 
ally increase its energy by increasing its line length or 
developing excitations, e.g. Kelvin waves. It is not clear 
which route to absorb energy would be favoured in 3D, 
although the orientation of the incoming sound is likely to 
play a deciding role. In systems of multiple vortices, the 
position of the vortices provides another route to modify 
the total vortical energy. Furthermore, symmetric exci- 
tations such as vortex rings and vort ex- ant i vortex pairs 
carry linear momentum and so can be expected to inter- 
act predominantly with sound waves of linear momenta, 
much like dark solitons. 
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VII. PROPOSED EXPERIMENTAL 
REALIZATION 

A. Set-up details 

The double trap used thus far (defined by Eq. is 
idealized and not experimentally realizable. Here we will 
demonstrate that the same qualitative phenomena occur 
in experimentally-realizable traps. We approximate the 
idealized double harmonic trap by considering a single 
elliptical harmonic trap which is split by a central Gaus- 
sian barrier. The "split trap" has the form, 



V(x,y) 



where, 



rn 
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(6) 



is an offset introduced such that the minimum of the 
potential is zero. The height of the barrier relative to 
the trap minimum Vo is related to these parameters via 
Vo = Vb — V m [ n . We set the barrier width to be d = 5£. 
Figure 10 compares the split trap of Eq. ([6| with the 
idealized double harmonic trap of Eq. ([2| . 

Experimentally, the harmonic trap can be formed via 
magnetic or optical fields, while the Gaussian barrier is 
formed by a blue-detuned laser beam aligned along the 
2/- axis. As previously, the system we simulate satisfies 
\i = lOhuo. Let us give some typical values for our units. 
We will assume some typical values: a 2D peak density of 
no = 10 14 m -2 , an axial trap frequency of uo z = 2tt x 1000 
Hz and a radial trap frequency of uj r = 2tt x 50 Hz. Then, 
for a 87 Rb ( 23 Na) condensate, the healing length becomes 
£ = 0.3 (0.7) /im, the speed of sound c = 3 (4) mms -1 
and the time unit (£/c) = 90 (175) /is. 



tion chosen. As such we will specify the ^-displacement 
of the vortices rather than the radial displacement. Fig- 
ure ITT] shows the phase diagram for one and two vortices 



in the split trap, as a function of the vortex position y v 
and the trap barrier height Vb- They show close agree- 
ment with the corresponding plots for the double har- 
monic trap, Figs. |2] and |6| a), and demonstrate the same 
regimes of the vortex dynamics, e.g. cross-over and in- 
ductive dynamics. 

Choosing a case from the two vortex system 
(Fig.JTTJb)) where the vortex dynamics is stable [y v = 3£ 
and Vo = 0.5/i ], we show the evolution of the vortices in 
Fig. [12] The vortices clearly drift in and out of the trap 
centre, out of phase with each other. This corresponds 
to the same transfer of energy observed in the idealized 
trap in Fig. [7Fa). The period of the energy transfer is 
around 1600(£/c) which is around 8 vortex precessions. 
Snapshots of the condensate density show the initial ap- 
pearance of the condensate (vortex 1 at large radius and 
vortex 2 at the trap centre), after a half-cycle of energy 
exchange (vortex 1 at the trap centre and vortex 2 at 
large radius) and after a full cycle of energy exchange 
(vortex 1 returns to large radius and vortex 2 returns to 
the trap centre). 



C. Precessing obstacle 

We now turn to the system with one vortex and a pre- 
cessing obstacle. The addition of a precessing Gaussian 
obstacle to a harmonic trap has been demonstrated ex- 
perimentally [9, 38 . The obstacle parameters we employ 
are obstacle width a = 2£, obstacle amplitude A Q b = /i, 
Vob — 3£, </>ob = and u; b = 0.234a;. The ensuing 

The vor- 



dynamics of the vortex is shown in Fig. 13 



tex gradually spirals into the trap centre, indicating the 



B. Cross-talk of two vortices 

Each well is no longer circularly symmetric and so the 
vortex displacement now becomes sensitive to the direc- 




FIG. 10: The experimentally realizable split-trap (grey solid 
line) and its vortex-free density profile (black solid line) along 
the y-axis. The analogous potential (grey dashed line) and 
density (black dashed line) for the double harmonic trap are 
shown. 



(a) 
10 

8 
6 
4 
2 




(b) 



0.5 



10 



+ + + + + + + 



1.5 



6 + 

+ 

4 + 
+ 

2 • 





0.5 



1.5 



FIG. 11: Phase diagrams for vortex dynamics in the split trap 
as a function of vortex displacement y v and barrier height Vb . 
(a) A single vortex in the right-hand well, displaced by y v 
from the well centre. Dots (crosses) correspond to the final 
state [after 5000 (£/c)] being a vortex (vortex-free) state, (b) 
A vortex in each well, of the same polarity. The final state is 
either two vortices (dots), one vortex (pluses) or no vortices 
(crosses). In (a) and (b) the light-shaded (dark-shaded) region 
corresponds to case I (II) of vortex dynamics. 
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FIG. 12: Vortex motion during cross-talk in the experimental 
split trap. Also showns are plots of the condensate density 
at times t = 0, 850 and 1550 (£/c). The vortices move anti- 
clockwise with respect to these density profiles. The color- 
density mapping is defined as black/ white=high/zero density. 
Parameters: Vb = 0.5/x and y v = 3£. 




FIG. 13: Vortex motion under parametric driving by a pre- 
cessing obstacle in the experimental split trap system. The 
condensate density is presented at times t = 0, 12000 and 
20000 (£/c). In the density plots, the hole in the left-hand 
well corresponds to the obstacle and the hole in the right-hand 
well to the vortex, and both move ant i- clockwise. Parameters: 
y v (t = 0) = 3£, Vb - 0.5//, a - 2£, A oh = /x, y oh = 3£ and 
cjob = 0.234cj). 



steady increase in its energy and angular momentum. As 
in the double harmonic trap, the vortex-obstacle energy 
transfer is slower than the vortex- vortex case, with min- 
imal position/maximal energy reached after ~ 10 4 (£/c). 
Note that the split trap system demonstrates the same 
resonant behaviour as in Fig. [9| for example, there is 
a sharp resonance in the obstacle precession frequency 
which closely coincides with the precession frequency of 
the unperturbed vortex (~ 0.234cj). 



D. Experimental observation 

The vortex-vortex and vortex-obstacle energy trans- 
fer could be observed through the evolution of the vor- 
tex position as it spirals in and out of the trap. Real- 
time tracking of vortices was recently demonstrated [13] . 
In brief this involves transferring (via pulsed microwave 
radiation) a small proportion of the BEC into an un- 
trapped state, allowing this representative cloud to ex- 
pand such that the vortex cores become optically resolv- 
able, and performing absorption imaging of the cloud. 
By repeating this at time intervals, the vortex motion 
can be tracked. For 87 Rb ( 23 Na) the period of the en- 
ergy exchange for the vortex- vortex case presented here 
is ~ 0.15 (0.3)s, while for the vortex-obstacle case it is 
~ 1.8 (3.5)s. These timescales are well within the ex- 
perimental lifetime of vortices in BECs, which can be as 
large as 10 s [7]. 

VIII. CONCLUSION 

We have shown that two vortices in a double trap 
can undergo a coherent cross-talk, periodically exchang- 



ing energy and angular momentum over long-range via 
sound waves. These observations are strongly analogous 
to the interaction of dark solitons mediated by linear- 
momentum-carrying sound waves. This adds further evi- 
dence for the striking acoustic similarities of vortices and 
dark solitons, despite no mathematically rigorous link. 
Sound waves artificially generated by a moving obstacle 
are similarly found to drive energy into a vortex. Cru- 
cially, for the sound to be absorbed by the vortex, it 
must carry angular momentum of the same sign as the 
vortex; if the sound has opposite angular momentum it 
is essentially invisible to the vortex. 

Our observations are robust: we have confirmed that 
we observe the same qualitative dynamics for different 
condensate parameters (e.g. a large, more Thomas-Fermi 
condensate) and double trap parameters. Importantly, 
our observations occur in experimentally realizable dou- 
ble well geometries, and within both the timescales of 
current vortex experiments and the timescale of other 
dissipation mechanisms, e.g. thermal dissipation [43| [44] . 
Our prediction of parametric driving of a vortex using a 
precessing obstacle could be used to counteract thermal 
dissipation of vortices, as suggested for dark solitons [18]. 

Trapping appears to promote sound absorption, sug- 
gesting that this may play a much greater role in trapped 
BECs than in homogeneous superfluids. An example 
may lie in the formation of vortex lattices in rotating 
atomic BECs. The nucleated vortices are observed ex- 
perimentally to crystallize into a vortex lattice, a pro- 
cess that requires dissipation. While this is provided by 
thermal dissipation at raised temperatures [45H47 : , ex- 
perimental and theoretical observations at very low/zero 
temperature suggest that crystallization is temperature- 
independent [48-51 . Vortex-sound interactions may pro- 
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vide such a mechanism. 

Establishing the rudimentary interactions between vor- 
tices and sound in quantum fluids will aid in exploring 
and understanding more complex scenarios, such as those 
present in large-scale quantum turbulence [3] and analogs 
of black hole superradiance based on scattering of sound 



waves from superfluid vortices [52j [53] . In future work 
we hope to explore sound absorption in 3D vortex lines 
and pursue acoustic ray models of sound absorption in 
trapped condensates. 

We acknowledge funding from the EPSRC. 
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